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»^ ■ Abstract 

A general construction principle of inflation rules for decagonal quasiperiodic tilings is 
^ . proposed. The prototiles are confined to be polygons with unit edges. An inflation rule 

OO . for a tiling is the combination of an expansion and a division of the tiles, where the 
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expanded tiles can be divided arbitrarily as far as the set of prototiles is maintained. A 
certain kind of point decoration processes turns out to be useful for the identification of 
possible division rules. The method is capable of generating a broad range of decagonal 



^ . tilings, many of which are chiral and have atomic surfaces with fractal boundaries. 

Two new families of decagonal tilings are presented; one is quarternary and the other 
ternary. Properties of the ternary tilings with rhombic, pentagonal, and hexagonal 
prototiles are investigated in detail. 



1. Introduction 

Since the discovery of quasicrystals (QCs) (Shechtman et ai, 1984), there has been 
an up-surge of interest in quasiperiodic tilings with non-crystallographic point symme- 
tries. Archetypes of such tilings in the plane are the Ammann-Beenker tiling (Griinbaum 
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&: Shephard, 1987; Beenker, 1982), the Penrose tilings (Griinbaum & Shephard, 1987; 
Penrose, 1974,1978,1979), and the Stampfli square-triangle tilings (Stampfli, 1986; 
Baake et al, 1992), which belong to the octagonal, decagonal, and dodecagonal Bra- 
vais classes, respectively. Their vertices form a quasiperiodic point set called a quasi- 
lattice (QL), which mimics the spatial arrangement of clusters in a two-dimensional 
(2D) QC. In real QCs, still a far wider variety of cluster arrangements are possible 
under different compositions and temperatures (Edagawa et al, 1994; Edagawa et al., 
2000). In order to advance our understanding of the real structures, it is therefore 
necessary to extend our scope for tiling models of QCs and to describe the structures 
in a more systematic fashion. In the present report, we mean by a 2D tiling a disjoint 
covering of the plane by edge sharing copies of a finite number of polygonal prototiles. 

In general, a 2D QL is constructed as a section of a 4D hypercrystal along the plane 
(physical space), where the hypercrystal is constructed from 'atomic surfaces' arranged 
periodically according to the relevant Bravais hyperlattice; each atomic surface extends 
only along the perpendicular space, which is an orthogonal complement to the physical 
space. The global characteristics or, more precisely, the local-isomorphism class of the 
QL is coded as the size and the shape of the atomic surface(s) (de Bruijn, 1981; 
Mackay, 1982; Kramer & Neri, 1984; Duneau & Katz, 1985). ^ 

In constructing 2D QLs by the section method, we would practically be restrained to 
polygonal atomic surfaces. However, certain important structures cannot be obtained 
in this way; the dodecagonal square-triangle tilings are typical examples which have 
rather complicated 'fractal' atomic surfaces (Stampfli, 1986; Baake et al, 1992; Smith, 
1993; Cockayne, 1994). Further examples with fractal atomic surfaces can also be found 
in the literature (Zobetz, 1992; Godreche et al., 1993; Niizeki, 2007a). These structures 

have been discovered by inflation methods, in which the shapes of atomic surfaces are 

Local isomorphism classes of QLs are further grouped into mutual-local-derivability classes (Baake 
et al, 1991; Baake & Schlottmann, 1997). 
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not given a priori. It is important to note that even for a known atomic surface 
with a fractal boundary, the computation of the structure by the section method is 
impractical because minute numerical errors cannot be avoided in judging which side 
of the boundary lies a given point near the boundary. Therefore, the section method 
does not suit the generation of such a QL, whereas an inflation method can work fine. 

So far, a systematic attempt has been made for generating ID binary tilings by 
inflation rules (Luck et al, 1993), revealing that the structures tend to have fractal 
(Cantor-set like) atomic surfaces. In more than two dimensions, however, inflation rules 
for quasiperiodic tilings have been rarely known except for those described above. This 
situation is simply caused by a hurdle in finding an inflation rule that does not produce 
any inconsistency throughout the entire structure. There is nevertheless some hope 
at present; a systematic inflation method for generating QLs, i.e. quasiperiodic point 
sets, in general dimensions has been developed recently (Niizeki, 2008). Let us call the 
latter method the point inflation scheme (PIS) and each of its inflation algorithms a 
point inflation rule (PIR). At this point, however, it should be remembered that a QL 
does not necessarily provide the vertices of a tiling. 

An advantage of employing tiling models for QCs lies in that they enable the entire 
structure to be decomposed into a finite number of prototiles, or fundamental struc- 
tural units. Another possible advantage may lie in their stronger geometrical con- 
straints, which are based on physical reasonings like the avoidance of unrealistic short 
distances. The present aim is therefore to develop a new inflation scheme for gener- 
ating quasiperiodic tilings. As in the case of archetypal tilings in the plane, inflation 
rules should proceed in the following steps: 1) an expansion of a tiling and 2) a division 
of the expanded tiles into tiles of the original size. Bear in mind that the expanded 
tiles can be divided in an arbitrarily way as far as the set of prototiles is maintained. 

In the present report, it is proposed that the division of the tiles is generalized with 
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the help of a point decoration process. The point decoration might simply be defined 
by a PIR, but this cannot prevent excess points from being generated in general. One 
needs to remove a part of the resulting point set so that the remaining points would 
be the vertex set of a tiling. Wherefore, the point decoration can be defined instead by 
taking an appropriate subset of the set generated by a PIR. In other words, the PIR 
is used to obtain the candidate positions for the decoration, where only a part of it is 
adopted. The present scheme turns out to be useful for discovering unknown tilings. 
For the sake of a self-contained presentation, the following arguments are confined to 
the decagonal case. Furthermore, the edges of any tiling to be considered are all given 
by one of the basic unit vectors, ej {j = 0, 9), which will be defined in Sec. 2. 

Sec. 2 is devoted to mathematical preliminaries, in which the decagonal Br avals 
module Aio is introduced; it is a projection of the 4D decagonal lattice onto the 2D 
physical space. The new scheme for generating decagonal tilings is presented in Sec. 3. 
In Sec.4, the present scheme is exercised and two new families of decagonal tilings 
are obtained; one is quarternary and the other ternary. The atomic surfaces of the 
new tilings are presented in Sec. 5. In particular, for the ternary tilings, the atomic 
surfaces with fractal boundaries are derived geometrically as the fixed sets of the dual 
maps associated with the relevant inflation rules. A statistical analysis is performed 
for the ternary tilings in Sec.6. After presenting further remarks in Sec.7, conclusions 
are given in the final section. 

2. Decagonal Bravais module 

Let us define the quasicrystallographic axes for a decagonal QL by the ten unit vectors 
Bj = (cos (jO), sin (jO)) with j = 0, 1, 2, 9 and 9 = it/5, pointing at the vertices of a 
regular decagon centered on the origin. Among the ten vectors, only four are linearly 
independent with respect to integral coefl&cients. A conventional set of basis vectors is 
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introduced as ej := e2j (j = 0, 1, 2, 3), pointing at four vertices of a pentagon from the 
origin, where the fifth vertex 64 := es is left unused because it can be represented by 
the four bases; 64 = — J2^=o ej. The basis set generates a Z-module of rank four called 
the decagonal Bravais module, denoted by the symbol Aiq. Remember that the point 
group of Aio is the dihedral group Diq. Furthermore, Aiq has an important property 
of scaling invariance rAio = Aio, where r = (1 + \/5)/2 is the golden mean (Niizeki, 
1989a). 2 

The decagonal Bravais hyperlattice £10 is defined in a 4D Euclidean hyperspace E4, 
and it is generated by the primitive basis vectors (ej,^ej-) {j = 0,1,2, and 3), where 
ej- := e^j (j^q^j ;^q) (Niizeki, 1989a; Yamamoto, 1996) and ^ is an arbitrary scale factor 
satisfying < ^ 7^ 1. The point group G of £10 is isomorphic to the dihedral group 
DiQ. The 2D physical space eH is an irreducible subspace of E4 with respect to G, and 
it is inclined against Ciq in an incommensurate fashion. The orthogonal projection of 
£10 onto E'I gives nothing but the decagonal Bravais module Aiq. 

The 2D perpendicular space E"*- is defined as the orthogonal complement to eH in 
E4; that is, E4 = E'I ® E-*-. Then the hyperlattice Ciq can be projected onto E"*- as 
well, generating another Z-module, K^q. Since the orthogonal projections from Ciq to 
both Aio and K^q are bijections, one can introduce a natural bijection tt between the 
Z-modules; Aj^q = ttAiq. 

We confine ourselves to the case when the vertex set Er of a tiling 7" is a subset of 
the Bravais module Aiq. In the ordinary cut-and-projection method (de Bruijn, 1981; 
Mackay, 1982; Kramer & Neri, 1984; Duneau & Katz, 1985), Er is generated as the 
orthogonal projection of (E'I + W) n £10, a cut of Ciq within a strip along E^, onto 
eII. The cross section of the strip, W(c E-*-), is equivalent to the atomic surface in the 

section method. The image of Eg- in E"*- is denoted as := txTit- In general, S;^ is 

^ In general, the ratio t is called the Pisot unit in the algebraic theory of Bravais modules. For the 
octagonal or the dodecagonal cases, it takes 1 + \/2 or 2 + ^3, respectively. 
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a dense subset of the atomic surface W, which on the other hand should be included 
in the closure of S^, i.e. C W C eJ. » 

3. Generalized point processes for tilings 

In the PIS (Niizeki, 2008) , a PIR is formulated as a set map that acts on an arbitrary 
subset L of Aio. proceeds by expanding L by a certain ratio a{> 1) and subsequently 
placing a copy of a certain motif S{c Aio) centered on every vertex. The ratio u can 
take any natural power of the Pisot unit r (the golden mean) associated with Aio, so 
that (jL remains to be a subset of Aiq. The motif 5 is a bounded set with a finite 
number of points comprised of one or more shells, each of which is an orbit of a point 
with respect to the point group Diq. One can put the procedure in a simple form as 
(t){L) := aL + S, in which the + symbol implies A + B = {a + b\^a e A,^be B}. Note 
that the resultant QL would have the point group Diq. 

The atomic surface of a QL generated by the PIS is identified as the unique attractor 
(Hutchinson, 1981) of the dual set map (j)-^ that acts in the perpendicular space; 
<p^{X) = a*X + (S-*-, where a* is the algebraic conjugate of a and = itS. (p-^ 
is nothing but an iterated function system (IFS), which is a common technique to 
generate fractal objects (Falconer, 1990). Indeed, the atomic surfaces of many QLs 
generated by the PIS have fractal boundaries. Various planar QLs have been found 
by the PIS (Niizeki, 2008) , which has also been applied to the case of icosahedral QLs 
(Fujita & Niizeki, 2008). 

Certain QLs can derive a tiling if the points are connected by uncrossed edges of 
a unit length, where the tiles are the regions bounded by these edges. Let us call 

this property the unit connectivity (UC) of the point sets. If a point set obtained by 

^ Some authors (Baake et al, 1991; Niizeki, 2008) prefer that an atomic surface W is a compact (i.e., 
closed and bounded) set in E""", hence, W = S:^. In this paper, however, more general language is 
sought so that a part of the boundary dW can be absent from W while fulfilling the glueing condition 
as explained in the text for the case of the RPH tilings, see Fig.8. 
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successively applying a PIR has the UC property, the PIR can be translated into an 
inflation rule for the tiling; an expansion of the tiles followed by their division into the 
original tiles. The division rules of the tiles are determined by the generated points 
in the tiles. In general, however, undesirable short distances prevent the generated 
point sets to be translated into tilings. Such redundancies can be discarded only by 
introducing an elimination step at each iteration to take an appropriate subset which 
fulfills the UC property. Then the expanded tiles are divided properly. We call the 
new inflation process a generalized point process (GPP). In the GPP, the role of the 
PIR is to generate the candidate positions for the vertices of a tiling. 

For now, the UC property is required by the point set at every iteration of a GPP. A 
GPP for decagonal tilings is formally given by the following steps: Step (I) an expan- 
sion of the tiling by the ratio a, where cr is a natural power of r (the golden mean), 
and Step (II) the decoration of every expanded tile by points, where the positions 
of the points should originate from Aio and are assumed to be determined uniquely 
by the shape of the expanded tile as well as those of its adjacent neighbours. In 
practice, candidate positions for the points of decoration are generated by a PIR, 
then an appropriate subset is taken to determine the division of the tiles. Note that 
there are degrees of freedom for the choice of the point decoration provided that it is 
determined uniquely within the first adjacent neighbors and that the set of prototiles 
is maintained. It is this generality that allows us to produce a broad class of new 
quasiperiodic tilings. The degrees of freedom can be promoted further by allowing 
the decoration to depend on farther neighbors. However, this would complicate the 
algorithm and will not be considered. 
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4. Examples of decagonal tilings 

Several decagonal tilings are generated by the GPP scheme. All the tilings presented 
below obey the UC property, in which the vertices are connected through the ten basic 
unit vectors Bj ( j = — 9) forming the edges of the tiles. Any vertex can therefore 
be written as I = noeo + niCi + 71262 + ^363, or alternatively it can be indexed as 
I = [nonin2n3]. 

4-1- Para-Penrose tiling 

The present tiling has been already reported (Niizeki, 2008) as one of the simplest 
decagonal tilings that can be generated by the PIS. The scaling ratio of the relevant 
PIR is cr = r^, while the motif S consists of two shells ([0000]) and ([1100]), where 
(l) := {gl\g G -Dio} is the orbit of the point I with respect to Diq. In addition to the 
four prototiles, i.e., the 36° rhombus (R), the regular pentagon (P), the crown (C), and 
the pentacle star (S), of the pentagonal Penrose tiling (PI) (Penrose, 1974; Griinbaum 
Sz Shephard, 1987), it has an additional prototile, namely the barrel shaped hexagon 
(H). This tiling is called the para-Penrose tiling by Niizeki (2008). 

In a single iteration, these tiles are expanded by the ratio and then a copy of S 
is placed on every vertex; see Fig. 1(a). One finds that each expanded tile is uniquely 
decorated with the points of 5's within its border. As the new points arc connected 
with unit edges, the expanded tile is divided into tiles of the original size, whereas in the 
vicinity of the boundary the tiles can be shared with the neighboring expanded tiles. 
In particular, when dividing an expanded P- or H-tile, segments near the boundary 
can be attributed to different kinds of tiles. This means that the division rule is not 
unique within the expanded tile itself Still, it turns out to be unique within the first 
adjacent neighbors as one can readily check in Fig. 1(b). Therefore, the present PIR 
satisfies the requisites for a GPP. 
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4.2. RPHC tilings 

The PIR for the para-Penrose tihng can be shghtly altered by removing a generated 
vertex inside every expanded C-tile at the symmetrical position. Obviously this does 
not affect the uniqueness of the decoration within every expanded prototile, and the 
resulting point set is shown to satisfy the UC property. As a consequence, the combi- 
nation of an S-tile and a P-tile (S-P complex) lying in the center of every expanded 
C-tile turns into one of a C-tile and an H-tile (C-H complex). Since the S-prototile 
may no longer appear except in the center of an expanded S-tile, it is a marginal one 
in the present case; that is, an iteration of the above GPP will generate a quaternary 
tiling with R-, P-, H-, and C-prototiles as shown in Fig.2(a). The division rule of every 
expanded tile depends on its first adjacent neighbors as in the case of the para-Penrose 
tiling. 

There are yet different ways to remove a vertex in an expanded C-tile without 
affecting the set of prototiles. Let us consider a gear-shaped complex (or a G-complex) , 
including an S-tile in the center and five adjacent P-tiles (see Fig.3(a)), lying on the 
bottom of the expanded C-tile. One or two of its five internal vertices can be removed 
so that the G-complex is to be divided in different manners. By removing one vertex, 
the G-complex is divided into one H-, one C-, and four P-tiles (Fig.3(b), left). There 
are five possible orientations, one of which has been taken in the preceding paragraph. 
Two other orientations are used in Fig.2(b) and (c), while the remaining two are their 
mirror images. Removing two internal vertices offers yet different ways to divide the 
G-complex. Of the five possible orientations, the only one without breaking the mirror 
symmetry is used in Fig.2(d). Two other choices are shown in Fig.2(e) and (f), while 
the remaining two are their mirror images. By removing two vertices, the G-complex 
is divided into one R-, two H-, and three P-tiles (Fig.3(b), right). 

We have constructed in total ten GPPs for tilings with the identical set of four 
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prototiles. Recall that the eight GPPs without mirror symmetry would generate chiral 
tilings. If the ten GPPs are applied in a mixed and arbitrary order, infinitely many 
tilings can be generated. Furthermore, if different GPPs are allowed to be applied at 
different locations at the same time, the possibilities become illimitable. 

4.3. RPH tilings 

Let us argue how the C-prototile can be excluded from the prototiles of RPHC 
tilings. Consider a turban-shaped complex (or a T-complex) composed of a C-tile 
and two adjacent P-tiles in an RPHC tiling. A T-complex is associated to every 36° 
angle of the expanded R- and C-tiles. The T-complex can be decomposed into three 
tiles, R-, P-, and H-tiles, by removing one of the two inner vertices, thereby breaking 
the mirror symmetry; see Fig.3(c). If all the T-complexes are (re-)divided in such a 
way, an RPHC-tiling turns into a tiling which does not include any C-tile. The two 
different ways to divide the T-complex are mirror images of one another, where left- 
and right-handed chiralities are associated with them. 

It is tempting to apply either one of the chiral rules for dividing all the T-complexes; 
then the two GPPs for both chiralities arc obtained. These GPPs maintain the set of 
the prototiles, R, P, and H. The left-handed GPP is shown in Fig.4, the repetition 
of which generates a tiling shown in Fig.5(a). The mirror images of these figures 
corresponds to the right-handed counterparts. Note that the point decorations within 
the expanded prototiles are not uniquely determined in this case, but they are unique 
within the first adjacent shell. The decoration uniquely determines the division of the 
expanded prototiles. The two GPPs may be applied in an arbitrary order, resulting 
in an infinite number of RPH tilings. A tiling shown in Fig.5(b) is an outcome of an 
alternate repetition of the right- and the left-handed GPPs, where the final GPP is the 
left-handed one. In both of the tilings shown in Fig. 5, a spiral pattern associated with 



lUCr macros version 2.1/31: 2007/06/16 



11 

the left-handed chirahty can be easily recognized. Again, much greater possibilities 
exist if the two GPPs are allowed to be applied at different locations at the same 
time. 

5. Atomic surfaces 

The atomic surface of a quasiperiodic tiling T is closely connected with the global 
characteristics of the structure. It can be inferred by projecting the vertices of a 
patch of T containing sufficient number of vertices onto the perpendicular space, E-*-. 
This is done in Fig. 6 for the six RPHC tilings and in Fig. 7 for the two RPH tilings, 
respectively. The convex hulls of all these atomic surfaces are the regular decagon 
whose vertices are given by ^Bj {j = 0, ...,9), corresponding to the atomic surface 
of the Penrose PI tiling. However, different types of erosion are observed near the 
boundaries even for tilings with the same set of prototiles. The symmetry of an atomic 
surface reflects that of the relevant tiling; all the atomic surfaces except the two with 
the mirror symmetric GPPs (Fig. 6(a) and (d)) have the lower point symmetry with 
the cyclic group Ciq. For each case, the erosion reveals a complicated fractal patterns. 
In particular, the three atomic surfaces that are shown in Fig.6(a),(c), and (e) exhibit 
hierarchical holes in the boundary regions. The rest of the atomic surfaces maintain 
the disc-like topology. 

Recall that, for a QL generated by the PIS, the atomic surface is simply identified as 
the unique attractor of the dual set map (j)'^{X) = a*X + S-^ in E-*-, where a* denotes 
the algebraic conjugate of a and X any subset of E"*- (Niizeki, 2008). For a tiling 
generated by the GPP scheme, on the other hand, the dual set map is associated with 
the determination of the candidate positions, from which a certain subset should be 
eliminated. The elimination is represented in E"*- as a subtraction of unnecessary parts 
from (f)^{X). In the following, the latter process is presented for the RPH tilings, in 
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which case geometrical rules to determine the atomic surfaces are identified. However, 
since such geometrical rules are not easy to generalize, the case of the RPHC tilings 
will be left open. 

The PIR for the para-Penrose tiling given in Subsec.4.1 gives S-^ = {([0000]), ([1010])} 
and a* = for the dual set map (j)-^. The fixed set of 0-*-, which is a moth-eaten 
version of the regular decagon (Niizeki, 2008), is the corresponding atomic surface. 
The small part that is left out from the decagonal atomic surface corresponds to a 
portion of concave vertices of C-tiles and S-tiles, resulting in the emergence of H-tiles 
in the para-Penrose tiling. 

In order to obtain the atomic surfaces for the RPH tilings, a subtraction process 
needs to be combined with the dual set map at each iteration. In Fig.8, the initial 
polygon Xq is defined as a decagonal star, which is known to be the atomic surface of a 
non-chiral RPH tiling (Papadopolos & Kasner, 2003). It is transformed by the dual set 
map to Xq = ^-^{Xq). The subtraction process for the chiral GPPs (Subsec.4.3) can 
be understood as a carving process near the boundary of Xq. Note that five strips are 
superposed on the next figure, while five more strips can be superposed up-side-down 
but are suppressed. These strips are cut by the boundary of Xq, while the two ends of 
each strip do not coincide through a translation. Accordingly, within each cut of the 
strips, two points can lie on a line parallel to the strip with the distance ^r. This leads 
to an excessive appearance of short distance 1/r in EH and hence to the existence of 
C-tiles. It turns out that for an RPH-tiling no such pair of points is allowed to exist 
within a single strip. Therefore, either or both ends of each cut must be carved so that 
they coincide through ^r-translations. 

The carving process for the left- or right-handed GPP is simply to carve a single end 
only. As shown in Fig.8, five strips are arranged to form a pentacle star, which can be 
traced either clock- wise or anti-clock-wise. Let us fix that the pentacle should be traced 
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clock-wise. Then for the left-handed GPP, the first end of each strip encountered while 
tracing the pentacle should be carved so that it coincides with the translate of the 
second end. These two ends fulfill the gluing condition; namely, if a point on one end is 
taken into account, the corresponding point on the other end should be discarded. The 
carving process for all the ten strips results in a new polygon Xi , which is the atomic 
surface of another RPH-tiling (chiral in this case). The carving process is denoted as 
7; , while the right-handed counterpart is denoted as 7^ in which case the opposite end 
of each strip is to be carved. 

The atomic surface of an RPH-tiling generated by the GPP scheme can be obtained 
by repeatedly applying 7/ ■ (p^ and/or 7,. • 0^. In the example shown in Fig. 8, only the 
left-handed process ji ■ (p-^ is used. In general, the subtraction process 7 associated 
with the removal of unnecessary points from the candidate positions can be rather 
complicated. At present, the author has been able to identify 7 only for the chiral 
RPH-tilings. 

6. Statistics of the RPH tilings 

The inflation matrices are introduced in the following to analyze the statistics of the 
RPH tilings. For each chirality, there are six types of tiles that are divided differently, 
as can be seen in Fig.4. In order to define the inflation matrix, the generated tiles need 
to be classified further into six types, so that their numbers can be counted in each 
expanded tile. Note, however, that the latter task depends on which chirality is to be 
used in the next iteration. In Fig.4, some of the small tiles are assigned different letters 
according to the chirality of the next GPP. There are four pairs of chiralities for the 
two successive iterations, //, rl, Ir, and rr. It is sufficient to consider two particular 
cases, II and rl, since the other cases are mere mirror images of the former two; the 
latter s have the same inflation matrices as the formers. 
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In the case U, the division rules shown in Fig.4 reveals that the six expanded tiles 
are divided as 

VA = Va + ^Vd + Fe, 

Vb = \Va + Vfe + 5y„ 

Vc = + 2Vfe + + (l - 27^) 
Vd = \Va + 3Vfe + + (3 - 2^^) Fd, 
F^, = 2y„ + 414 + (4 - j ^d, 
Vf = Wa + 2Vb + Wc + 2Vd, 

where Vx represents the volume of the tile labeled 'x'. In the case Ir, the corresponding 

formulae are 

VA = Va + ^Vs+Vf, 
Vb = \V-a + H + 5Fe, 

Vc = 2^^ + Vfc + 4y, + (1 - ^) v-d, 
Vd = \V-a + 2H + 3Fc + (3 - 2^^) Fj, 
Fij = 2l^s + 2V'5 + 2Fg + (4 - % 
= 3^^ + 2H + 4Fe + 2Fj. 

It follows that the inflation matrices for all the four cases are given by 
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One can check that the maximal eigenvalues of both Mi and M2 are r^, correspond- 



ing to the rate of volume increase under a single GPP iteration. The relevant right- 
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eigenvectors are common to the two matrices, reading 




= (0.236,0.691,0.691,1,1,1)*, (5) 

where the superscript t indicates the transposition. These six components provide the 
relative area (volume) of the tiles labeled 'a' to '/'• 

On the other hand, the left-eigenvectors for the common maximal eigenvalue give 
the number ratio of the tiles. For Mi and M2, the relevant left-eigenvectors are 

ui = (r^r3 + l,2T^r3- 1,1,0) 

= (6.854,5.236,8.472,3,236,1,0), (6) 

U2 = (r^r3-l,2T^ + 2,r3-l,0,l) 

= (6.854,3.236,10.472,3,236,0,1), (7) 

respectively. That the eigenvectors are different might seem troublesome if the two 
GPPs are applied in an arbitrary order. It turns out, however, that no problem is 
caused by the difference, since the inflation matrices multiplied by the wrong eigen- 
vectors will just give the right ones, i.e., 

W1M2 = T^U2, (8) 
U2M1 = T^Ul. (9) 

It follows that the number ratios of the tiles depend only on the chiralities of the final 
two iterations. The above two distributions for the six types of tiles do not cause a 
difference in the number ratios of the three prototiles, since the sum over components 
for P-tile or H-tile are common to the two left-eigenvectors. It follows that the number 
ratios of the three prototiles in the present RPH-tilings are R : P : H = 1 : 2 : 1/t, 
while the mean volume of the tiles is l/r times that of an H-tile. 
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The differences in the left-eigenvectors manifest themselves in the statistics of local 
arrangements of tiles in the relevant tilings. For instance, the second components of 
ui and U2 representing the relative frequencies of P-tiles labeled '6' are different. This 
is manifested in the frequencies of the local centers of five-fold symmetry in the two 
tilings shown in Fig. 5; the difference can be rather significant from the viewpoint of 
the structural stability as well as the physical properties if these tilings are to be 
used for modeling physical QCs. Importantly, the difference should also be connected 
to the boundary shapes of the atomic surfaces, which point will be left for a future 
investigation though. ^ There are on the other hand three types of local centers of two- 
fold symmetry (Niizcki, 1989b) located (i) at the centers of R-tilcs labeled 'a' or 'a', 
(ii) at the centers of H-tiles labeled 'e', 'e', or '/'j and (iii) at the mid-edge positions 
between two adjacent P-tiles both labeled '6'. The frequencies of these two-fold centers 
are common to the two cases. 

7. Further remarks 

The generality of the GPP scheme is capable of generating a number of unknown 
quasiperiodic tilings, many of which has a chirality. Let us briefly consider how the 
structure factor of a chiral tiling looks like. In our examples, the breaking of the mirror 
symmetry is carried by a small part of the atomic surface near the boundary, while the 
main body of the atomic surface maintains the mirror symmetry. Since the latter part 
provides the main contribution to the structure factor, the chirality is only manifested 
in relatively weak Bragg peaks. A structure factor is shown in Fig.9 for the tiling 
shown in Fig.5(a), assuming a point scatterer on every vertex. 

In the GPP scheme, the removal of unnecessary points from the given candidates 

"* Identifying the classes of local centers of symmetry (Niizeki, 1989b) in quasiperiodic tilings is an 
important task, which has been addressed in a systematic way for the cases with polygonal atomic 
surfaces (Niizeki, 2007b). A generalized argument is hence necessary to handle the cases with fractal 
atomic surfaces. 
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for the point decoration forms a critical step. This makes the GPP the most general 
and robust technique for generating decagonal tilings. Furthermore, the application 
of the basic idea to the octagonal as well as the dodecagonal cases is straightforward. 
Indeed, for most of the known tilings, whether the atomic surfaces are polygonal or 
fractal, inflation rules can be re-phrased as GPPs, i.e. the combination of an expansion 
step and a point decoration step. 

Let us consider the particular case of the dodecagonal square-triangle tilings. In this 
case, each of the vertices in the expanded tiling is decorated by a three shell motif S, 
comprising the origin, an inner hexagon, and an outer dodecagon. The hexagon can 
take two different orientations, which can be chosen at random (Smith, 1993; Gahler, 
1988) or according to a deterministic rule; for instance, the local coordination can be 
used to fix the orientation (Hermisson et al, 1997). One can see a similarity of the 
situation to the case of decagonal tilings that are presented in this paper. 

Different tilings that are generated by applying different GPPs in different orders 
are all members of the random tiling ensemble with the same set of prototiles. They 
are likely to have energies very close to each other, so do they have similar statistical 
weights. Since they form a deterministic subset of the relevant random tiling ensemble, 
the relevant contribution to the entropy is called the 'deterministic entropy', a term 
coined by A. P. Smith (Smith, 1993). 

Either of the left-handed and the right-handed GPPs for the RPH tilings is defined 
by exclusively applying one of the two division rules for the T-complex; see Fig.3(c). 
One readily recognizes that these two division rules are mutually connected through a 
phason flip involving three tiles. Accordingly, the RPH tilings embrace an abundance 
of flipping sites. This has a signiflcant implication when a physical realization of a 
similar structure is to be considered, since the phason degrees of freedom must play 
an important role in the structural stabilization. One should also bear in mind that 
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the same kind of phason flips has been observed in situ in a d-Al-Cu-Co QC at 1123 
K with a transmission electron microscope (Edagawa et al., 2000). 

8. Conclusions 

A general inflation scheme for generating decagonal quasiperiodic tihngs has been 
proposed. In the new scheme, inflation rules are practicafly comprised of three steps: 
an expansion with the ratio a, a decoration of every vertex by a finite motif S, and an 
elimination of unnecessary points by local rules. At every iteration, the resulting point 
set should be unit-connective and form the vertex set of a tiling. Since the point decora- 
tion process can be readily generalized, various new tilings are expected to be found. 
The usefulness of the present scheme has been demonstrated by generating several 
new decagonal tilings, among which the family of ternary tilings has been analyzed 
in detail. The concept can be exported not only to the octagonal and dodecagonal 
cases, but also to the icosahedral cases (including P-, F-, and I-type Bravais classes), 
for which only a few tilings have been known. 

The author is greatly indebted to K. Niizeki for instructions in the mathematics 
of quasilattices and for a comment on the issue of local centers of symmetries in the 
present tilings. He is also grateful to T. Ogawa for a comment on spiral patterns in 
the chiral RPH tilings. 
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Fig. 1. (a) The point decorations of the expanded prototiles for the para-Penrose tihng 
are determined by superposing decagons centered at ah the vertices, (b) A square 
patch of the para-Penrose tihng with expanded tiles being indicated with the thicker 
hues shows the division rules for the expanded tiles. 
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Fig. 2. Square patches of RPHC tilings generated by GGPs with expanded tiles being 
indicated with the thicker lines. They are distinguished by the division of the G- 
complex lyings on the bottom of every expanded C-tile. The division rules of the 
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expanded P- and H-tilcs arc also affected by adjacent C-tiles. (a-c) A single point 
inside every G-complex is eliminated, while the mirror simmetry is retained only in 
(a). 
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Fig. 3. (a) One G- and three T-complexes associated with an expanded C-tile (the 
hatched area) are shown, (b) Depending on whether one or two internal vertices 
are removed, a G-complex is divided in two different ways, which may further take 
five distinct orientations relatively to the expanded C-tile. (c) A T-complex can be 
divided in two ways depending on the choice of an internal vertex to be removed. 
In (b) and (c), the black dots indicate the removed vertices. 
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ABC 




D E F 

Fig. 4. The left-handed GPP for the RPH tihngs is shown. An expanded P- or H- 

tilc is divided in several different ways depending on its first adjacent neighbors. 
Expanded tiles are thus classified into six classes according to their division rules, 
and they are labeled 'A' to 'F'. The black dots indicate the removed vertices inside 
T-complexes. The generated tiles can also be classified into six classes but in two 
different ways depending on the chirality of the subsequent GPP (see also Sec. 6). 
For instance, a P-tile labeled '(6, c)' implies that it is classified into the 6-subclass 
if the next GPP is left-handed while into the c-subclass if the next GPP is right- 
handed, where the bar indicates the mirror inversion. A mutually inverted pair is 
denoted by a starred label; that is, c* = (c, c) and d* = {d,d). The label 'a*' (i.e. 
'(a,a)') for each R-tile is suppressed merely for the reason of space. Tiles with the 
label '/' only appear if the two successive GPPs have opposite chiralities, while tiles 
with the label 'e' (or 'e') only appear if the same GPP is applied successively. 
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Fig. 5. Square patches of RPH tilings generated by repeating the left-handed GPP 
(a) and by alternating between the right- and left-handed GPPs (b). Spirals in the 
arrangement of P-tiles are emphasized with the gray spots. The black dots indicate 
the centers of the expanded P-tiles with the five-fold symmetric division rule (Fig.4- 
B), showing that the frequencies of local five-fold centers are somewhat different 
between the two. 
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(a) . - , (d) 




Fig. 6. (a)~(f): The atomic surfaces for the RPHC tihngs in Fig.2(a)~(f), respectively. 

The gray areas represent the projections of a large patch containing over 250,000 
vertices onto E"*-. Fine details including hierarchical holes (pits) are visible. 
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Fig. 7. (a) and (b): The atomic surfaces for tlic RPH tilings in Fig. 5(a) and (b), 
respectively. The gray areas represent the projections of a large patch containing 
over 240,000 vertices onto E-*-, while the boundary lines are obtained with the 
geometrical rules as presented in Fig.8. 
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Fig. 8. The dual map for the left handed GPP is illustrated. Starting from a star 
decagon Xq for a non-chiral RPH tiling, the dual set map c/)-^ and the subtraction 
process 7; is alternately applied. For applying the dual set map, a copy of the 
reduced figure, e.g. a*Xo, is placed on every point of indicated by the black 
dots. For 7;, a single end of all the ten strips are carved as described in the text. 
Only five strips are depicted by the thin lines. By iteration, the composite of the 
two maps, 7; • (p^, generates a series of figures Xj {i = 0,1,2,...) (thicker solid 
lines), each of which corresponds to a ternary RPH tiling. The limit figure X^o 
is the boundary of the atomic surface W for the RPH tiling generated solely by 
the left-handed GPP. Note that a triangular region close to the boundary of W is 
divided into a central regular pentagon (A), 3/5 of which is occupied, and the two 
adjacent isosceles triangles (B), half of which in total is occupied. 
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Fig. 9. The structure factor of the RPH-tihng generated by the left-handed GPP 
depicted in Fig. 4. The area of each spot corresponds to the intensity. The chirahty 
is manifested in the weakest spots. 



Synopsis 

A general construction principle of inflation rules for decagonal quasiperiodic tilings is pro- 
posed. 
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